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Abstract 

The Dicke spin-boson model is composed by a single bosonic mode and an ensemble of TV 
identical two-level atoms. Assuming thermal equilibrium with a reservoir at temperature 
we consider the situation where the coupling between the bosonic mode and the atoms generates 
resonant and non-resonant processes. The thermodynamic of the model is investigated. Next 
we introduce dipole-dipole interaction between the atoms. We investigate the transition from 
fluorescent to super-radiant phase and the quantum phase transition in a situation where the 
dipole-dipole interaction between the atoms generates entangled states in the atomic system. 
We proved that, the critical behavior is not modified by the introduction of the dipole-dipole 
interaction. 
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1 Introduction 



According to classical physics, the knowledge of the state of a composed system yields complete 
knowledge of the individual state of the parts. This is the principle of separability [1] [2]. For 
example, if we known the state of a classical system of particles, we known the state of each 
particle. However, this is not true in the context of quantum mechanics. One of the consequences 
of the superposition principle is the introduction of the concept of entanglement. After interacting, 
two quantum systems can end up in a non-separable state, i.e., states that can not be factorized 
into a product of the states of its sub-systems. These states of such composite system are called 
entangled states [3]. For instance, in a bipartite entangled state, each part losses its quantum 
identity. Entanglement leads to the non-locality properties of quantum mechanics, and action at a 
distance at a speed greater than the speed of light [4] [5]. Discussing the hidden variable theories 
based on local realism, Bell has shown that quantum theory is supposed to be non-local [6] [7] 
[8]. It was be shown that any pure entangled state of two spin-^ violates corresponding Bell's 
inequality [9] [10] [11] [12]. 

With the development of quantum information [13] [14] [15] and its application in computation 
and communication [16] [17] [18] [19] [20] [21] [22] [23], the entangled states have been attracted 
enormous interest. For instance, new tests of quantum mechanics can be implemented using some 
entangled states, as for example the Bell and the GHZ states [24] [25]. Also, several quantum 
protocols can be realized exclusively with the help of entangled states [26] [27] , as for example the 
quantum dense coding protocol proposed by Bennett and Wiesner [28], or the quantum telepor- 
tation protocol proposed by Bennett and co-workers [29]. In the teleportation of entanglement a 
single qubit in an arbitrary state can be transferred, exchanging information over long distances, 
without physically transferring the system itself. The maximally entangled pair of qubits is cru- 
cial in the quantum teleportation. Actually, there are many physical systems that can be used 
to implement a quantum computer with the quantum logic gates. Some experimental devices are 
based on cavity quantum electrodynamics, trapped ions and nuclear magnetic resonances. For 
example, for the realization of quantum logic gates, a many-body system prepared in an entangled 
state has been proposed [30]. Another promising system used to implement quantum computation 
is the quantum dot array proposed by Loss and DiVincenzo [31]. The basic problem that arises 
in this area of research is how to create entangled in many-body systems [32] [33] [34] and also 
how to generate systems which are not affected by the environment, overcoming the problem of 
decoherence. 

In the simplest case of two atoms, it is quite important to demonstrate creation of entanglement 
on such system. There are in the literature different methods for detection of entangled states 
of two interacting atoms. One method is based in measure the angular intensity distribution of 
the fluorescent field emitted by these two atoms, since it is well known that radiation emitted by 
atoms exhibits directional properties [35] [36] [37] [38] [39] [40]. The other is based in the dynamic 
of the population inversion of the system. For instance, the properties of spontaneous emission 
from two identical entangled atoms interacting with the modes of a bosonic field was investigated 
by Guo and Yang [41]. They have shown that the time evolution of the population inversion, 
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which is proportional to the radiation intensity, depends on the degree of entanglement of the 
initial state of the system. 

Suppose that in a many-body system an entanglement has been created in a given portion of 
such system. A fundamental question is how to known whether a state of a many-body system is 
entangled or not [42]. Before continue, we would like to briefly describe the theoretical machinery 
that are needed to define entanglement in a pure or mixed state. To define separability and 
entanglement in bipartite systems, we use the properties of the state vector or density operator. 
Since a pure state is entangled if it is not separable, it can be shown that a pure state is separable 
if and only if the reduced density operators of sub-systems represent pure states. Therefore a 
pure state is entangled if and only if the reduced density operators for the sub-systems describe 
mixed states. In practice, a pure state is separable if and only if the quantum fluctuations of all 
observable (linear self-adjoint operators acting in the respective Hilbert space of each sub-system) 
are uncorrelated. If at least one pair of linear self-adjoint operators have correlated quantum 
fluctuations, the pure state is entangled. For mixed states we define a separable mixed state if 
the two sub-systems have the same purity and von Neumann entropy. An entangled mixed state 
is one that is not separable. We would like to stress that experimental tests for separability and 
entanglement of pure states using correlated quantum fluctuations are not available for mixed 
states. Strong correlation between two operators does not means that there is entanglement for a 
system described by a mixed state. 

Critical phenomena is a cooperative effect characterized by fluctuations in the order param- 
eter, where all scales of length are important. For instance, the two-point correlation function 
(the expectation value of products of local observables) near the phase transition does not de- 
crease exponentially, but polynomially, as the distance between the points goes to infinity. The 
system exhibit long range order. Since, for pure states the correlated quantum fluctuations of 
two observables defines an entangled state, we should expect a close connection between second 
order phase transition at zero temperature and entanglement. We would like to note that recently 
Emary and Brandes [43] [44] discussed the connection between quantum phase transition and the 
chaotic behavior that emerges in the full Dicke model for finite N, where the energy level-spacing 
statistics changes from Poissonian to one described by Gaussian ensembles of the random matrix 
theory [45] [46]. 

Let us now focus our attention for the super-radiant phase transition. Consider an ensemble 
of iV two-level atoms, all prepared into the excited state. Each atom can emit a photon by 
spontaneous emission. In the situation where there is not coupling between the atoms, each 
atom radiates independently. This is called the fluorescent phase. Let us consider the following 
experimental environment: the iV two-level atoms are in an optical cavity where all atoms are 
in resonance with a single mode of the field. For an ensemble of atoms in a volume with linear 
dimensions small compared to the emission wavelength, they start to radiate spontaneously faster 
and strongly than the ordinary fluorescent phase. In this situation the radiation rates becomes 
quadratic dependent on the number of atoms. This cooperative process is called super-radiance. 
Recently Lambert and co-workers [47] investigated the entanglement properties of an ensemble 
of N two-level atoms interacting with a field mode in the super-radiant phase [48] via the von 
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Neumann entropy. The maximum entanglement occurs near the critical region. The authors 
also studied the system at zero temperature and how does the entanglement is affected by a 
quantum phase transition that occurs in the system [49] [50]. By the other hand, in a system 
of iV atoms, the presence of the dipole-dipole interaction can generates entangled states in the 
many-body system. See for example Ref. [51], where was discussed in detail how to generate multi- 
particles entanglement of atoms trapped by harmonic potential, interacting with a classical field 
and between them via the dipole-dipole interaction. Therefore one can imagine that if entangled 
has been created in a portion of the many-body system, the critical behavior of the system changes 
due to the strong correlation between the atoms. 

The aim of this work is to investigate if the entanglement between two-level atoms, generated by 
the dipole-dipole interaction, is able to change the critical properties of the system. For instance, 
it is well known that the dipole-dipole decohering effect is able to suppress the super-radiant 
emission [52] [53]. Note that we are not able to quantify the degree of entanglement present in 
the model since we are using functional integration methods to investigate the thermodynamic of 
the model. First, we investigate the full Dicke model where the quantum phase transition and a 
phase transition from fluorescence to super-radiant phase, at some temperature (3~ l , in the system 
of N atoms interacting with a bosonic field is analyzed. Second, introducing the dipole-dipole 
interaction we discuss the effects of entanglement in this many-body system, showing that the 
critical temperature of the transition from the fluorescent to super-radiant phase is not modified 
in this case. We proved that the spectrum of the collective bosonic excitations of the model is 
unaffected by the dipole-dipole interaction. In this case, there is also a quantum phase transition, 
at some values of the physical parameters of the model. For the reader interested in the the study 
of entanglement in many-body systems close to the quantum phase transition, see for example the 
Refs. [54] [55] [56]. 

Recently, using the path integral approach with functional integration method the analytic 
behavior of thermodynamic quantities in the full Dicke model [57] [58] [59] [60] [61] [62] was 
presented [63] . See also the Refs. [64] [65] . The full Dicke model is similar to the Jaynes-Cummings 
model [66] where the non-resonant processes in which the atom and the field are excited or de- 
excited simultaneously, known in the literature as the anti- Jaynes-Cummings model. In Ref. [63] 
the study of the nonanalytic behavior of thermodynamic quantities in the full Dicke model, allow 
the authors to evaluate the critical transition temperature. It was shown that the system develop 
a quantum phase transition and a phase transition from ordinary fluorescence to super-radiant 
phase at some critical temperature. Based in this analysis, the spectrum of the collective bosonic 
excitations, for different situations were presented. It was study the rotating-wave approximation, 
the counter rotating-wave approximation and also the general case. As we discussed, the rotating- 
wave approximation ignores energy non-conserving terms in which the emission (absorption) of a 
quantum of a quantized field is accompanied by the transition of one atom from its lower (upper) 
to its upper (lower) state. In a situation where only non-resonant processes contribute, the full 
Dicke model present a second order phase transition from the ordinary fluorescent to the super- 
radiant phase respectively, at some critical temperature and also a quantum phase transition, 
i.e., a phase transition at zero temperature. In the last case there are no thermal excitations, 
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therefore the phase transition is driven by the quantum fluctuations. The interesting result is the 
fact that it is possible to have a condensate with super-radiance in a system of N two-level atoms 
coupled with one mode of a bosonic field where only non-resonant processes contribute. 

In the Ref. [67], still neglecting the direct interaction between the atoms, it was considered 
also two different models, assuming that a single quantized mode of a bosonic field interacts with 
a ensemble of N identical two-level atoms. Again we assume that the system is in thermal equi- 
librium with a reservoir at temperature Analytic properties of the partition functions of the 
models were also investigated. First, it was study a modified version of the model discussed by 
Chang and Chakravarty, Legget and others [68] [69] [70] [71], which has been used to analyze 
dissipation in quantum computers. Owing to the coupling between the two-level systems and a 
bosonic reservoir, this model presents destruction of quantum coherence without decay of popu- 
lation. Since the interaction Hamiltonian of the model generates a quite particular non-resonant 
processes, it is shown that the partition function is analytic for all temperatures, and therefore 
there is no second order phase transition in the model. Second, it was investigated a model where 
the coupling between the bosonic mode and N two-level atoms is intensity dependent [72] [73] 
[74] , introducing also the couplings g\ and g 2 for rotating and counter- rotating terms respectively. 
At low temperatures, the contribution coming from the counter-rotating terms dominates over the 
rotating ones. When the coupling constant is given by g 2 = (ujq Q) 2 , where Uq is the energy of 
the single-mode bosonic field and Q is the energy gap of the atoms, a quantum phase transition 
appears. 

In this paper the physical system in consideration is an ensemble of two-level atoms interacting 
with a single mode of a bosonic field. It is a natural question to ask if non-Gaussian terms can 
change the critical behavior of this spin-boson model. We study two models, the full Dicke 
model and the same model where we include the dipole-dipole interaction term between the two- 
level atoms. We proved that the critical temperature that characterize the phase transition from 
fluorescent to super-radiant phase is the same in both models. Also the quantum phase transitions, 
at some values of the physical parameters of both models are the same. We are using the path 
integral approach with the functional integration method to investigate the thermodynamic of the 
models, which is given by the analytic properties of the partition function. The paper is organized 
as follows. In section II we present the fermionic full Dicke model and the model with the dipole- 
dipole interaction. In section III we discuss the functional integral for the full Dicke model. In 
section IV we study the full Dicke model with the dipole-dipole interaction. Conclusion are given 
in section V. In the Appendix the theory of pure and mixed ensembles and the reduced density 
operator is briefly presented. In the paper we use ks — c — h — 1. 
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2 The N two-level atoms-Bose field interaction Hamilto- 



nians 



In this section we consider a very general situation where the system under investigation 
contains a large number of two-level atoms. In order to describe the dynamics of the the two-level 
atoms and the bosonic mode we have to introduce the Hamiltonian governing the interaction of 
the quantized Bose field with free atoms [75] [76] . Free means that there is no interaction between 
the atoms. Therefore let us consider a Bose quantum system B, with Hilbert space HS B ^ which 
is coupled with iV atoms, with Hilbert space TiS^ . Let us assume that the whole system is in 
thermal equilibrium at temperature The Bose quantum system is a sub-system of the total 
system living in the tensor product space ® TC^\ Let us denote by Hb the Hamiltonian 

of the quantized Bose field, by Hq the free Hamiltonian of the N two-level atoms and Hi the 
Hamiltonian describing the interaction between the quantized Bose field and the N two-level 
atoms. The Hamiltonian for the total system can be written as 



where Ib and Iq denotes the identities in the Hilbert spaces of the quantized Bose field and the 
N atoms. In the Eq. (1), 5 is a small coupling constant between the atoms and the quantized 
Bose field. 

The main purpose of this section is to discuss iV-atoms-Bose field interaction Hamiltonian. 
Therefore, let us introduce the Dicke operators to describe each atom. The free j — th atom 
Hamiltonian will be denoted by Hjp , since we are using the Dicke representation. Therefore, we 
have 



where | i )j are orthogonal energy eigenstates accessible to the j—th atom and ui\ are the respective 
eigenfrequencies. Using Eq. (2) and the orthonormality of the energy eigenstates we can write the 

j — th atom Hamiltonian H$ as 



H = H B <g> Iq + Ib ® H Q + 8H h 



(1) 




(2) 



2 




H^\\i)(i\)r 



(3) 



i=l 




(4) 



+ 



(|2)(1|) 



(5) 
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(6) 



(7) 



and finally 

%)= (|l)(2|) r 

Combining Eq. (3) and Eq. (4), the j — th atom Hamiltonian can be written as 

ti D - W a {j) + -I u l +uj 2 I, 
where the energy gap between the energy eigenstates of the j — th atom is given by 

n^ = u^-u^. (8) 

Shifting the zero of energy to + u;^) for each atom, the j — th atomic Hamiltonian given 

by Eq. (7) can be rewritten as 



(9) 



Note that the pseudo-spin operators <r^, and <t^ satisfy the standard angular momentum 
commutation relations corresponding to spin | operators, i.e., 



and finally 



00' 



-a 



or 



(10) 

(ii) 

(12) 



The model that we are interested is composed by the Hamiltonian of the atoms, with the contri- 
bution of the quantized Bose field Hs, and the interaction Hamiltonian The Hamiltonian of 
the total system is given by 



H = I B <g) H Q + H B <g> I Q + SHj = 

oo oo N oo 

i B ® n E <?tj) + Y,wkblh ® /q + ^=EE ^ ( & * + & I) ® (^5) + %■)) • (1 3 ) 



j=i fe=i v -/V J=1 fc=1 

where the first summation in the right hand side is 



N 



E ^ 4) = °!d 



i + ... + 1 ® i 



1® fiW o-fjv), 



(14) 
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and 1 denotes the identity in the Hilbert space of each two-level system. The quantity g^j is the 
coupling coefficient. Each two- level system is described by a Hilbert space 7i^ x \TiS 2 \ HS N \ The 
composite atomic system is described by the iV-fold tensor product space Tip) ® <S> H^... 

One point which is important to stress is that in the Eq. (13) there are terms so called counter- 
rotating terms. In the rot at ing-wave- approximation we ignore energy non-conserving terms in 
which the emission (absorption) of a quantum of a quantized field is accompanied by the transition 
of one atom from its lower (upper) to its upper (lower) state. As we discussed before, the rotating- 
wave-approximation ignores terms in which the j—th atomic raising (lowering) operators multiplies 
the field creation (annihilation) operator. 

We can also introduce an atom-atom interaction, which is relevant in the study of entangled 
states. In an entangled system, the state of the composite system can not be factorized in to 
a product of the states of its sub-systems. For example in the case of two two-level atoms, the 
interaction Hamiltonian of two atoms including the dipole-dipole interaction reads 

1 2 

H (ii) = m E H (H) a t) ®%y ( 1 5 ) 

Let us note that the interaction Hamiltonian given by Eq. (15) is a particular case of the quantum 
Hamiltonian describing an axially symmetric ferromagnet, with the coupling constants J?- and . 
In terms of the pseudo spin operators we have 



H = ~z2 U °h °h + -h + %)°u) J • ( 16 ) 

<i,j> \ v 7 / 

The < i,j > symbol denotes a summation over pair of first neighbors. In the absence of the 
"dipole-dipole" interaction the pure Hilbert space of the two atomic system is spanned by the 
states | g± ) <S> \ <? 2 ), | gi ) <8> | e 2 ), | e± ) (g) \g 2 ) and | e x ) <g> | e 2 ), where g and e denotes respectively 
the ground and the excited state of the two atomic system. If we include the dipole-dipole 
interaction term in the form of Eq. (15), the vectors \g±) ® | e 2 ) and \e±) <E> | gi ) are not more 
eigenstates of the Hamiltonian of the atomic system. It can be shown that these two vectors states 
must be substituted by the two entangled states, known in the literature as maximally entangled 
states [12] [38] 

l s > = ^7| d e i) ® U2> + Ui> ® |e 2 » (17) 

and 

|a) = "y| (|ei> ® \g 2 ) ~ \gi) ® |e 2 ». (18) 

The Hilbert space of the system is spanned by the collective states of the two identical interacting 
two-level atoms given respectively by (\gi) <S> \ gi), \ s), \ a), \ ei) <S> |e 2 )). As stated before, 
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both the states | s ) and | a ) tell us nothing about the state of the first or second atom separately. 
Therefore when we include the dipole-dipole interaction between the atoms, we open the possibility 
of realize entangled states. 

In the spin-boson model with the dipole-dipole interaction included, we are considering the 
question of how does the dipole-dipole interaction change the critical temperature, where again 
the system exhibit a phase transition from fluorescence to super-radiance. The Hamiltonian for 
the spin-boson model with the dipole-dipole interaction reads 

1 N 

^3 

N oo N oo 

I B n 4) + E b t b k ® J Q + -J^ T E E 9ik (b k + b{) ® (<r+ + a 0) ) . (19) 
j=l k=l v iv j=1 k=1 

Each two-level atom interact with all other atoms of the ensemble with the same coupling strength, 
i.e., if(jj) = A, and the summation is over all the atoms. This is a model with one infinite range 
interaction. This model is based in the models describing spin systems with long range interaction 
[77] [78]. In the next section we discuss first the full Dicke model, showing the presence of a 
quantum phase transition and also a super-radiant phase at some critical temperature. 



3 The functional integral for the fermionic full Dicke model 

The typical situation in cavity quantum electrodynamics is the case where a single atom is 
coupled to a cavity mode. Now suppose N identical two-level atoms in an optical cavity whose 
linear dimensions are small compared to the wavelength associated to the bosonic mode. This 
ensemble of identical N two-level atoms interacting linearly with one mode of a bosonic field 
defines the full Dicke model. The Hamiltonian of the system reads 

JV n a N 

H = Is ® E 9 + & f b ® I B + -j= £ (b + &t) g, (,+ + . (20) 
j= i z ViV j=1 

In the above equation g is the coupling constant between the atoms and the single mode of the 
bosonic field. The b and W are the boson annihilation and creation operators of mode excitations 
that satisfy the usual commutation relation rules. 

The aim of this section is to prove that a model with an interaction Hamiltonian where we 
take into account non-resonant processes in which the two-level systems and the bosonic modes 
are excited or de-excited simultaneously, presents a phase transition from normal to super-radiant 
state at some temperature with the presence of a condensate and also a quantum phase transition 
at some critical coupling. We should mention that the results of this section are not original, but 



8 



by completeness we include in the paper. The Euclidean action of the model is S, where Hp is 
the full Hamiltonian for the fermionic full Dicke model. This quantity can be written as 



it / \ 
H F = ^b*{r)b{r) + - ^(a*(r)a t (r) - ^(r)^(r)j + 

+ j= E(aKr)p i (r)b(r) + a< (r) P*(t) b*(r) ) + 

+ E («i(r))9;(r)6(r) + a?(r) ^(r) 6*(r)). 



(21) 



Note we are introducing two coupling constants, #i and f° r the rotating and counter- rotating 
terms, respectively. As we discussed before, the main reason for this is that we are interested 
in to identify the contribution of the real and virtual processes in the phase transition with the 
formation of the condensate. For this purpose we must to find the partition function of the system, 
therefore we calculate the formal quotient -J^, being defined by the equation: 



Z 



f[drj\e 



So 



(22) 



where the functional integral is with respect to the complex functions b*(r) and b(r) and Grass- 
mann Fermi fields a*(r), a«(r), (3*(r) and A(r). Here the expression [drj] is the functional measure. 
It is possible to show that an expression for which is given by [63] [67] 



Z 



a(0) + 2c(0) 



a(0) - 2c(0) 



n 

u> >0 



1 - a(w)) - a(-uj)^ - Ac 2 (uj) 



■1/2 



+ 



+0(N- 1 ), 

where a(u) and c(u) in the above equation are given, respectively, by 



a{uj) 



and 



c(u) 



' g\ (tt - iu) _1 + gl(Vt + iu) _r 
(u - iu) 

( 9x92^1 



tanh 



pn 



tanh 



pn 



(23) 



(24) 



(25) 
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To obtain the thermodynamic limit we have to take the limit (N — > oo) in Eq. (23). We turn 
out to the discussion concerning the local elementary excitation of the ground state. To find the 
energy spectrum of the bosonic collective excitations we have to use the equation 



4c 2 (cu) 



a{uj) 



a{ —uj) 



0, 



(26) 



and making the analytic continuation {iuj — > E), we obtain the following equation 



1 = 



9t 



ilu; 



e 2 )(q 2 
( i 



E 2 ) 



{u 2 - E 2 ) 



+ 



tanh 2 
1 



4 



+ 



4ft' 



(ft — E) 2 (ft + E) 2 (ft 2 - E 2 ) 2 



tanh' 



+ 



g 2 (n - E)-i + g 2 (n + E)^ + g 2 (n + E ) _1 + g 2 (n - E ) _1 



(lu - E) 



(u> + E) 



+ 



tanh 



Solving the above equation for the case j3 1 = j3 c 1 we find the following roots 

E l = 

and 

E2 = 



g 1 (Q + ^ ) 2 + 92 (0 - ^o) : 



1/2 



(27) 
(28) 
(29) 



V (91 + 92) J 

Its low energy state of excitation is a Nambu-Goldstone mode, since the continuous U (1) symmetry 
has been broken. Therefore by using the non-relativistic version for the Goldstone theorem, the 
energy spectrum can not have a gap above the ground state. Now, let us present the critical 
temperature and the energy level spectrum of the collective bosonic excitations of the model with 
the rotating-wave approximation, where g\ 7^ and g 2 = 0. The result obtained by Popov and 
Fedotov [79] [80] is recovered, where the equation 



a(0) 



and 



tanh 



(¥)-'■ 



gives the inverse of the critical temperature, (5 C . It is given by 



P c = — arctanh( -^-] . 



\ 9i 



(30) 
(31) 

(32) 
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The order parameter of the transition is the expectation value of the number of excitation as- 
sociated to the bosonic mode per atom, i.e., lim N ^^^p- ^ 0. Note that again cuq, f2 and g\ 
define also a non-zero critical temperature where the partition function is no more analytic. We 
may expect a super-radiant phase for the temperature (3~ l given by Eq. (32). The energy level 
spectrum of the collective Bose excitations in this case is 

£i = 0, (33) 

and 

E 2 = tt + uj q . (34) 

In this case, there is also a quantum phase transition, i.e., a zero temperature phase transition 
when the coupling constant g\ attains the value g\ = (uj Q)^. Now we will show that is possible 
to have a condensate with super-radiance in a system of N two-level atoms coupled with one mode 
of a Bose field where only virtual processes contribute. In the pure counter-rotating wave case, 
i.e., g\ = and g 2 7^ 0, the inverse of the critical temperature, j3 c is given by 

P c = — arctanhf— ^- V (35) 
l I \ g 2 J 

and the spectrum of the collective Bose excitations given by 

£i = 0, (36) 

and 

E 2 = \Q-u \. (37) 

Again, there is a zero temperature phase transition when g 2 = (cuo Q) 2 . In this case there is no 
thermal excitations and the phase transition is driven by the quantum fluctuations. A comment is 
in order concerning the spectrum of the Bose excitations. In both of the cases: working with the 
pure counter-rotating or the rotating-wave terms, there is a phase transition. In the case of the 
rotating-wave approximation g± ^ and g 2 = 0, there is a Nambu-Goldstone mode (E — 0). In 
the pure counter-rotating case gi = and g 2 7^ also there is a Nambu-Goldstone (gapless) mode. 
We show that it is possible to have a condensate with super-radiance in a system of N two-level 
atoms coupled with one mode of a Bose field where only virtual processes contribute. Since the 
energy E 2 0, local elementary excitations of the ground state with low energy can easily be 
created causing a significant fluctuation effect. Unfortunately we are not able to evaluate these 
effect in the systems. 
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4 The functional integral for the spin-boson model with 



dipole-dipole coupling 

The aim of this section is to prove that the introduction of the dipole-dipole interaction does 
not modify the temperature of the phase transition from normal to super-radiant state and the 
presence of a condensate. Also, the energy spectrum of the collective bosonic excitations of the 
model is the same as in the full Dicke model. Introducing the dipole-dipole coupling in the full 
Dicke model, where an ensemble of identical N atoms interacts linearly with one mode of a bosonic 
field, in this section we use the same approach of functional integrals to study the thermodynamic 
of the model [63] [67] [79] [80]. We would like to stress that in this new situation, the Grassmann 
functional integrals are not Gaussian. In order to perform the integrals we introduce auxiliary 
variables in the calculation. The total Hamiltonian of this spin-boson model is defined by: 

1 N 

N N 

I B <g) Q ]T a z {j) + uj fct b <g> I Q + -j= J2 (b + 6 f ) ® (<tJ) + % } ) ■ (38) 
j=i V -/V j=1 

In the above equation the quantity g is the coupling constant between the atoms and the single 
mode of the bosonic field. The b and are the boson annihilation and creation operators of mode 
excitations that satisfy the usual commutation relation rules. As usual, we are using the pseudo- 
spin operators af-y a7-s and afo which satisfy the standard angular momentum commutation 

relations corresponding to spin \ operators. We are also shifting the zero of energy + w 2 ) for 
each atom and defining Q = cu 2 — 

Let us define the Fermi raising and lowering operators a], ati, (3\ and that satisfy the anti- 
commutator relations ctjO;] + ajctj = and fiifij + pjfii = Sij. We can also define the following 
bilinear combination of Fermi operators, ajctj — at\(3i and finally (3\ai. Note that <7^, at^ and 
(7 7j\ obey the same commutation relations as the above bilinear combination of Fermi operators. 
Therefore, we can change the pseudo-spin operators of this spin-boson model by using the bilinear 
combination of Grassmann Fermi fields 

a z (l) — (don - A f A) , (39) 
(7+ — . ajPi , (40) 



and finally 



(3}a t . (41) 
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Using the Eq. (38), the Euclidean fermionic action for the model can be written as 



v Ufm + £ l Ut) ^a + ff(T) %z | , _ L dTHp(T) . (42 ) 




The Euclidean action S is given by Eq. (42) where H F is the full Hamiltonian for the fermionic full 
Dicke model with dipole-dipole interaction. As in the previous section, introducing two coupling 
constants, g\ and g 2 , for the rotating and counter-rotating terms, the Hamiltonian becomes 

H F = u b*(T)b(r) + - £ («*(tWt) - ^(r)^(r)j + 

i N N , s 

+ ^E H (iJ) a *i(r)Pi(r)P*(T)a j (T) + ^ £ (a*(r) /3,(r) b(r) + a< (r) ^(r) 6*(r) ) + 



+ -4= E («i(r))9;(r)6(r) + a?(r) /3,(r) 6*(r) ) . (43) 



Let us define the formal quotient of two functional integrals, i.e., the partition function of the 
interacting model and the partition function of the free model. Therefore, we are interested in to 
calculate the quantity given by Eq. (22) where S = S(b,b* ,a,a^ , f3, ft) is the Euclidean action 
given by Eq. (42). Let us define So = So (6, b*, a, or ', /3, ft) which is the free Euclidean action for the 
free single bosonic mode and the free two-level atoms and finally [drj] is the path integral measure. 
Note that in Eq. (22) we have functional integrals with respect to the complex functions b*(r) 
and b(r) and Grassmann Fermi fields cx*(t), CKi(r), ft(r) and A(t). Since we are using thermal 
equilibrium boundary conditions in the imaginary time formalism, the integration variables in 
Eq. (22) obey periodic boundary conditions for the Bose field, i.e., b{j3) = 6(0), and anti-periodic 
boundary conditions for the Grassmann Fermi fields i.e., aii((3) = — aj(0) and Pi(/3) = — A(0). The 
free action of the bosonic field is given by 



S (b,b*) = J q dr(b*(r)-^--u; b*(r)b(r)^ . (44) 

In order to obtain the effective action of the bosonic mode we must integrate over the Grassmann 
Fermi fields in the expression given by Eq. (22). The problem now confronting us is the fact that 
in the action given by Eq. (43) there is a non-Gaussian term, corresponding to the dipole-dipole 
interaction. Although it is not possible to integrate this term directly, using auxiliary variables of 
integration we can circumvented this difficult. The following integral is helpful 



I N \ - N 

N 



^ x*A jA x, I = (detA)- 1 / f[^^e^'-^ A ^ + ^^^ + ^^=^^ ! , (45) 

■> i=i — 27ri 
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where this integral is valid for commutating variables. To continue we have to study the non- 
Gaussian contribution to the partition function. In terms of the Grassmann variables we define R 

as 



R = exp|-l ^pr tt ;(r)ft(r)ff^;(r)a 3 (r) 

= ex P E jf dr {PXtWt))* H tj /3;(r)a,-(r) j , (46) 

and since each pair of Grassmann variables and commutes, we have (aifl*)(aj(3*) = 
(ctj(3*)(ai(3*). Therefore in the functional version of the model we can substitute the Gaussian 
integral given in Eq. (45) in Eq. (46). We get 

R = a (detHr J[dr]e^ tf*-WV^>+^ Jo"* + , (47) 

where the auxiliary fields rj(r) satisfy periodic boundary conditions, i.e., rj(0) = rj(/3) and also 
r*(0) = r*(/3). In the above expression the quantity ao is a numerical factor that can be absorbed 
in the normalization factor and the term [dr] = Y\f [dr^\ [dr*] is the functional measure for the 
auxiliary complex fields. Substituting the Eq. (47) in the partition function Z, we obtain 

Z = a (det H)- 1 J [dr]e-££=i /o" dTr . T W fr « 1 ^« / [^] , (48) 

where in Eq. (48) we have S r = S r (ri,r*,b,b*,ai, ...). Consequently the last functional integral 
J[drj] e Sr depends on the auxiliary fields r*(r) and rj(r). Moreover, the action S r can be separated 
into a free action for the bosons and a Gaussian fermionic part. The action S r can be written in 
the form 

rf> N + 

S r = S (b, b*) + drJ2 Pi(r) M r Pi {r) . (49) 

We would like to stress that in Eq. (49), instead of writing M(r i: r*, b*, 6), for simplicity we use 
the notation M r , i.e., M r = M(r i: r*,b* : b). The column matrix Pi(r) is given in terms of the 
Grassmann Fermi fields 



P\(r) = (/3*(r) a*(r)) (50) 
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and the matrix M r is given by 



M r 



( 



\ 



-N- l / 2 ( gi b(r)+g 2 b* (r)-ri(r) 



d _ 



(51) 



/ 



The complex functions and Grassmann Fermi fields r^r), b(r), ol^t) and /3j(r) can be represented 
in terms of a Fourier expansion. Therefore, we have 



n(r) = /r 1/2 £^)^ T , 

V 

&(T)=/r 1/2 £&M^ T , 



and 



PiM=/r 1/2 £pi(p)e ipr . 



(52) 
(53) 

(54) 



Since the complex functions rj(r) and 6(r) obey periodic boundary conditions, and the Grassmann 
Fermi fields a^r) and /^(r) obey anti-periodic boundary conditions, we have that v = ^p, to = ^jp 

and p = ( 2n + 1 ) 7r ; for integer ra. They are respectively the bosonic and fermionic Matsubara 
frequencies. Substituting the Fourier expansions in the action given by Eq. (49) we get 



v 



S r = J2(iv- u ) b*(u) b(u) + £ £ Pi(p) m pqPM) , 

u) p,q j=l 



(55) 



where the matrix M pq = M pq (ri, r*, b, b*) is given by 



M, 



pq 



(ip + Q/2)5 pq -(N(3)-^Ql q 

-(Npy^Qpg {i P -n/2)s pg 



and the function Q pq = Q pq (ri, b, b*) can be written as 

Q P q = 9i Kp ~q)+92 b*(q - p) - n(p - q) . 



(56) 



(57) 



Since the integrals with respect to the Fermi fields are Gaussian, we may integrate over these 
Grassmann variables. This procedure yields 



r ( N \ N 

/ [ dri(p)} exp( £ £ Piip) M pc PM) = II det M * > 

\p,q i=l > j=l 



(58) 
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where the matrix Mj = Mj(rj, r*, 6*, 6) is a block matrix of the following form 

m,J ' P+ ^! -w)- 1 ^ 1 ). (59) 

In the above equation / is the identity matrix and the components of matrix P are P pq = pS pq and 
Q was defined in Eq. (57). Our aim is to present the quotient -J^, defined in Eq. (22). Combining 
the results obtained in Eq. (48), Eq. (49) and Eq. (58) we get that is given by 

(detH)- 1 j[dr]e^i=^ r > )H ^ rj[v) /[^(&)] e EJ"o) b *H b H n ^ det M . 
J[d V {b)}e^^-^ b '^ b ^det N M{0,0) 
where the functional measure [dr)(b)] is defined by 

[dr]{b)]=l[db{uj)db*(uj) . (61) 

Ul 

In the Eq. (60) the matrix M(0, 0) is given by 

In order to simplify the calculations, let us change variables in the following way: 

/ \ 1/2 

&H- ? 6( W ) (63) 

and 

/ \ 1/2 

&»- ? &». (64) 

Note that Eq. (64) is not the complex conjugate of Eq. (63). It is not difficult to see that after 
these changes of variables, the denominator of the Eq. (60) turns out to be equal to unity 

|[d77(fe)]exp(-7r^6*HfeH) =1. (65) 
We can express the ratio M- by the integral 



z 

J- = (dettf)- 1 J [dr] e^iE^W^^iW J [ drj ( b )] exp (s eff (b) ) , (66) 
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where after doing the Fourier transform for the auxiliary fields the functional measure becomes 
[dr] = nf^ EL dr,i(uj)dr*(uo). The effective action of the bosonic mode, S e ff(b) is given by 

N 

S eff = -7t ]T b*(uj)b(uj) + £ In det (J + A t ) . (67) 

u> i=l 

The matrix in the determinant of the above equation is given by 

det(I + Ai) = det^M- 1/2 (0,0)Mi(6*,6)M- 1/2 (0,0)y (68) 
Performing the product of matrices we obtain the matrix A afterwards. Therefore we can write 

The components of the matrices Bi and Cj are given by 

m Pq = -(^) 2 (^-|) 2 (^+|) 2 x 

x _ r * (g _ p) ) (70 ) 

\^u> -i(q-p) yju -i(p-q) ) 

and 

(QU = (^) 2 (^-f) 2 (^+§) 2 x 

X ( ^ 9 ' Hp - q) + V ; 9 ^ {q - P) -r^-q)). (71) 
V ^u -i(p-q) ^Juj -i{q-p) J 

In order to perform the functional integral given by Eq. (66) we must to find a manageable 
expression for det (I + A). For this case we can use the following identity 

det (/ + Ai) = det (/ + Bid) -> exp^r^C^)) . (72) 

So using this approximation given by Eq. (72) we can find an expression for the ratio defined 
in Eq. (66) of the following form 

— = (det//)-' A^^E.-W^-'C*) A)-JW §1, (73) 
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where 

^1 = f [d v (b)] e~ * ^ b * H ( X " aH ) 6(w) + * ( b(w) c{w) b{ - w) + b * (w) c(w) 6 * ( - w) ) 
Z J 

xe f E„ (di(«)6(w) + d2(«)6*(w)) 

In the Eq. (73) the coefficients a(u) and c(u) are given respectively by 

= ^(n-M-'+aKn + M-) tanh 

\ (uo - iu) J V 4 / 



X 



diH = ~ , ; 1 ( 7T^ E + E ^(-a;)! tanh f £5 ) 



and 



d 2 (w) = . f7T^-E^H + 7^— E r ?(-^)) tanhf^V 

Performing the Gaussian integral in the bosonic variables in the Eq. (74), we finally obtain 

— = / xexp — J2 Tl ; tt , ,x — E^Hr w + 

\iV 2 V (l -a(-u)) (l -a(w)J -4c 2 (cj) ^ 



+ E t: : * - 2 , : (E^H^(-) + E r *H r ] 

w ( 1 — a(— a;) J ( 1 — a (a;) J — 4c 2 (uj) 



where Jo, A\(u) and ^2(0;) are defined respectively by 

(l - a(0) - 2c(0)) 1/2 (l - a(0) + 2c(0)) 1/2 

1 

(l - a(-w)) (l - a(uj) - 4 c 2 H ' 



tanh f ^ 



Al(a) = (— gl— + — ^— - tanh (gi? ) - 2 - f . V 
7r(ft-iu/) V^o-^ u; + «u; V 4 / (o;g -o; 2 )(SZ + za;)y 

= 2WD2 i 2\ taIlh \~T J ■ 

7T (cJq — ^ )(ii + CJ 2 ) \ 4 / 
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Comparing the partition function for the full Dicke model in Eq. (23), with the partition function 
of the model with dipole-dipole coupling in Eqs. (73) and (79), we can see that the poles in 
the partition function are the same and therefore the temperature of the phase transition from 
normal to super-radiant state and the presence of a condensate is not affected by the dipole-dipole 
coupling. The energy level spectrum of the collective bosonic excitations of the model is the same 
as the original full Dicke model. The quantum phase transition is not modified in this situation. 
An important question here, is why the full Dicke model with the dipole-dipole interaction has 
the same critical behavior as the original model? 



5 Conclusions 



With the development of quantum information and its application in computation and com- 
munication the entangled states have been attracted enormous interest, since several quantum 
protocols can be realized exclusively with the help of entangled states. For example new features 
in cryptography becomes possible utilizing the non-classical properties of entangled pairs of par- 
ticles. The basic problem that arises in this area of research is how to measuring entanglement 
in many-body systems. The Schmidt decomposition can be used to measure correlation between 
two sub-systems in a joint pure state. An important property is that the reduced matrices (see 
appendix A) of both sub-systems written in the Schmidt basis are diagonal and have the same 
positive spectrum. For more than two entangled sub-systems the Schmidt decomposition is in 
general impossible to be implemented. For mixed states of two or more subsystems there is not a 
Schmidt decomposition, and different measures of entanglement are nonequivalent. The structure 
of entanglement in many-body systems is much more difficult to study then for bipartite systems. 

In a system of N atoms, the presence of the dipole-dipole interaction generates entangled 
states. If entangled has been created in a portion of the many-body system, the critical behavior 
of the system can changes revealing the strong correlation between the atoms. In this work we 
investigate if the entanglement between the atoms, generated by the dipole-dipole interaction, 
changes the critical properties of the system, as the critical temperature that characterize the 
phase transition from fluorescent to super-radiant phase and also the spectrum of the bosonic 
excitation of the model. 

First we consider the full Dicke model composed by a single bosonic mode and an ensemble of N 
identical two-level atoms. Assuming that the system is in thermal equilibrium with a reservoir at 
temperature we consider the situation where the coupling between the bosonic mode and the 
atoms generates resonant and non-resonant processes. Secondly we consider the full Dicke spin- 
boson model with the dipole-dipole interaction between the atoms. We show that, introducing the 
dipole-dipole interaction term, the critical temperature that characterize the phase transition from 
fluorescent to super-radiant phase does not change. Also, the spectrum of the bosonic excitations 
of the model is unaffected by the dipole-dipole interaction. Finally, the quantum phase transition 
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is not modified by the dipole-dipole interaction. Furthermore, as discussed previously there is a 
zero temperature phase transition when g 1 = and g 2 ^ ()■ In this case there are no thermal 
excitations and the phase transition is driven by the quantum fluctuations. An important point is 
that there is a close connection between quantum phase transition and entanglement. A quantum 
phase transition can be characterized by long range correlations. As the same way, the definition 
of entangled pure state is that in the system at least one pair of the observables has quantum 
correlations. 

We might ask why we have obtained this result, since we expect at principle that the non- 
Gaussian terms in the Hamiltonian can modify phase transition properties. Remember that the 
spin-boson Hamiltonian is a particular case of the Hamiltonian describing a system of bound 
charges and a radiation field in the dipole approximation [81], where we are assuming that the 
system is confined in a region whose linear dimensions are small compared with the wavelength 
of the bosonic mode. In the full Dicke model, field non-uniformity is disregarded. All atoms 
are in the same environment. The introduction of the dipole-dipole interaction with infinite 
range does not modify this situation. Each two-level atoms is coupled to infinitely many others. 
The above discussion leads us to conclude that we are solving mean-field models. Other heuristic 
argument is related to the fact that mean field theories give exactly the right description of critical 
properties for large d, where d is the dimensionality of space. Since our results are independent 
of the dimensionality, we are solving mean field models. That is the reason why the dipole-dipole 
interaction is not able to modify the critical properties for this cooperative phenomena of super- 
radiance. 

We conjecture that the critical properties are different if we go beyond the mean field approx- 
imation. First let us assume that linear dimensions of the system are not small compared with 
the wavelength associated to the bosonic mode. In this case we have to take into account the 
spatial variation in the coupling coefficient between the pseudo-spin operators at-y aZ\ and the 
boson annihilation and creation operators of mode excitations. We note that Li and co-workers 
[82] investigated the quantum phase transition and the super-radiant phase in the full Dicke model 
assuming that the dimensions of the atomic ensemble are much larger than the wavelength of the 
bosonic mode. This system has a critical temperature of the phase transition from normal to 
super-radiant state and also the presence of a condensate. With the introduction of the dipole- 
dipole term, the critical temperature of the super-radiant phase transition must change. This 
subject is under investigation by the authors. 
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A The theory of pure and mixed states and the reduced 
density operator. 



In the theory of quantum systems, the notion of composite quantum systems is fundamental. 
The Hilbert space of a composite quantum system is the tensor product space of the Hilbert spaces 
describing its sub-systems. Being more specific, let us consider a bi-partite quantum system S 
made of two quantum systems and S^, with respective Hilbert spaces TC^ and . The 
state space H of the combined system (S^ 1 ) + S^) is given by the tensor product of the Hilbert 
spaces pertaining to the sub-systems and S^ 2 \ Therefore we have that the total Hilbert space 
can be written as 

n = n {1) <g> n {2) . (a.i) 

If we take fixed orthonormal basis {| (p i ^ )} and {| ip ^ )} in TiS 1 ^ and respectively, a general 
state in the tensor product space TC may be written as 

l^> = E«*;k« (1) > ®kj (2) >- ( A -2) 

where {| (p i ^ ) ® | (p j ^ )} is a basis of the tensor product space. 

Suppose that we have an operator A W acting only in TiP^ and an operator A ^ acting only 
in H^K Therefore we define the tensor product A' 1 ' <g> A ^ showing how does this quantity acts 
on any element of the basis of the tensor product space. We have 

(A W ® A (| <p t W ) ® | Vj W )) = (A W I (1) )) ® (A ( 2 ) | ^ ( 2 ) )) . (A.3) 

This important definition lead us to its linear extension for an arbitrary state ■ Using the Eq. 
(A.2) and the Eq. (A.3) we get 

(A W ® A ( 2 )) | V» ) = £ a „ (A W I Vi W ) ) ® (A ( 2 ) | ^ W )) . (A.4) 

We also have that any product of operators acting in Ti can be represented as a linear combination 
of tensor products. We have 

A = Y, 0<* (1) ® A a ^y (A.5) 

a 

The observable of the sub-system S ^ take the form A' 1 ' ® I ( 2 \ while observable of sub-system 
take the form I W ® A( 2 ). 
Suppose that we are interested in the behavior of the sub-system S^ 1 ). This behavior is de- 
termined by the expectation values of the operators { }, where the label (1) means that each 
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operator acts on the states of the subsystem alone, and the index a = 1,2, ...,N means 
operators in . In other words, the label a denotes the family of commutating operators 
( [X a ^\X/3 = o) , acting in the state \ip). Since the state vector | ip ) of the combined system 
can be expanded using the orthonormal basis { | ip i ^ ) } and { | (p j ^ ) } we get for any operator 
{X^ }, is its expectation value given by 

(X«> = U|A (1) hA>. (A.6) 

Note that we drop the indices a in the operator. Using the Eq. (A. 2), the fact that the operator 
X* 1 ) does not act on the vectors that span the Hilbert space H^ 2 \ and the orthonormality of the 
{ | (p . ( 2 ) ) } basis, we get 

(^ (1) > = E«^-<^ (1) l^ (1) k 4 (1) >- (A.7) 

ijk 

Defining = J2j a ij a tji we have that the expectation value of the operator X^ in a general 
state \ip ), defined by ( X^ ), is given by 

(A (1) ) = E^(^ (1) I^ (1) I^ (1) ), (A.8) 

ik 

or 

(X«) =tr[X( 1 V (1) ], (A.9) 

where 

p (1) = E^ (1) X^ (1) I- (A.io) 

ik 

The quantity p^ is called the density operator of the sub-system The density operator 

describes the state of the system interacting with other system in the same way as the state vector 
describes the state of an isolated system. 

Is a fundamental interest to measuring entangled in many body systems, where we have usually 
mixed sates. For simplicity, let us consider a two-component quantum system S^ 1 ) and 
in a pure state A pure state is separable if and only if the reduced density operators 

represent pure states. Essentially separable states satisfy the classical separability principle. In a 
bipartite quantum system with an entangled state | <E» S ) , measurements in its sub-systems show 
that | $ s ) contain information in the measurements in S^ 1 -* and and also correlation between 
the measurements. To test whether the sub-systems and are entangled or not, we perform 
the Schmidt decomposition of the state vector | $ ). If the Schmidt number is greater then one, 
the state | $ ) is an entangled state. Therefore a pure state is entangled if and only of the reduced 
density operators for the sub-systems describe mixed states. Using the Schmidt decomposition it is 

possible to evaluate the density operators p^ ixed = 2V(2) I $ ) ( $ I and p^ ixed = 7V(i) | <&)(<& | . The 
two reduced density operators have the same eigenvalues, and the two sub-systems have identical 
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von Neumann entropies. On the other hand, if the Schmidt number is one, the reduced density 
operators for or represent pure states and the bipartite system is not an entangled state. 
The von Neumann entropy of a pure state is zero. Suppose the operators and defined 
as X « <g> I & an d / (i) <g> x ^. The fluctuation operator are defined as 6 = X« - ( X«) 
and 5 X^ = X^> — (X^). In practice, a pure state is separable if and only if the quantum 
fluctuations of all observable (linear self-adjoint operators acting in the respective Hilbert space 
of each sub-system) are uncorrelated. If at least one pair of linear self-adjoint operators have 
correlated quantum fluctuations, the pure state is entangled. 

Note that the von Neumann entropy attains the maximum S(p) = InN, in the case of com- 
pletely mixed state where p = jj X) 4 (ipil, where N is the dimension of the corresponding 
Hilbert space. Thus quantity is a monotonically function of the degree of entangled between a 
pair of sub-systems. Therefore the von Neumannn entropy of the density matrix for S^ 1 -* or 
of the bipartite system is a measurement of entangled [42]. 
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